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Abstract 

In this paper we construct some natural geometrical objects on the 1- 
jet space J 1 (R, R 3 ), like a nonlinear connection, a Cartan linear connection 
(together with its d-torsions and d-curvatures), a jet "electromagnetic" d- 
field and its geometric "electromagnetic" Yang-Mills energy, starting from 
a given dynamical system governing the three-dimensional motion of the 
knee in the mathematical model introduced by Grood and Suntay. The 
corresponding Yang-Mills energetic surfaces of constant level (produced 
by this knee dynamical system) are studied. 
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1 Short introduction 

The 1-jet spaces are basic objects in the study of classical and quantum field the- 
ories (see Olver [10]). For such a reason, a lot of authors (Asanov [2], Saunders 
[T2] and many others) studied the differential geometry of 1-jet spaces. Using as 
a pattern the Miron-Anastasiei's Lagrangian geometrical ideas (see [5]), Balan 
and Neagu have recently developed the single-time Lagrange geometry on 1-jet 
spaces (see [3]), which is very suitable for the geometrical study of the solutions 
of a given ODE system, via the least squares variational method initiated by 
Udri§te (see [T3]). 

In the present days, the biomechanics of the lower limb is an extremely 
fruitful field of research. This research field is of interest for a lot of biomech- 
anists (e.g., see Barbu [4] and Ro§ca [lljl. In this direction, it is important to 
note that the differential geometrical methods from 1-jet spaces are very use- 
ful for studying biomechanics (see Ivancevic [7]). In such a perspective, this 
paper applies the jet geometrical results from monograph [3J to the Euler's dy- 
namical equations that govern the three-dimensional motion of the tibia with 
respect to the femur, in the biomechanical and mathematical model introduced 
by Grood and Suntay [5], and studied further by many scholars (e.g., see Hcfzy 
and Abdel- Rahman [B]). 



2 Jet single-time Lagrange geometry produced 
by a non-linear dynamical system 

Let us present now the main jet single-time Lagrangian geometrical results that, 
in our opinion, may characterize a given non-linear dynamical system. In this 
way, let us consider the jet fibre bundle of order one J X (R, R") ->Kx R ra , where 
n > 1, whose local coordinates (t, a; 4 , x\), i = 1, n, obey the rules 



t = t(t), x i =x i (x j ), x\ 



dx l dt 
dxj~dl 



Let X = (xi l X(x k )j be an arbitrary distinguished (d-) tensor field on the 



1-jet space J 1 (K, R n ), whose local components transform by the rules 



x (i) = 



dx l dt 



X 



The d-tensor field X produces the jet first order ODE system (Jet dynamical 
system) 

x\ = X®(x k (t)), yi=T~^, (I) 

where c(t) = (x l (t)) is an unknown curve on R" (i.e., a jet field line of the 
d-tensor field X), and we have 



dx l 
~dt' 



V i = l.n. 



Let us consider now the jet least squares Lagrangian function (attached to 
the dynamical system ([1]) and to Euclidian manifolds (R, 1) and (R n ,5y)) 



JLS 



ODEs 



J 1 ! 



expressed by 



JLS 



(x k ,Xi) — Sij x\ X ^(x) x{ X ^(x 



EH 



X^(x) 



(2) 



where x — (x k ) k= —- It is obvious that the global minimum points of the jet 
least squares energy action 



.-ODEs,- 



(c(t))= 3LS uubjS (x k (t),x k (t))dt, te[a,b], 



2 



are exactly the solutions of class C 2 of the jet dynamical system (p}. In other 
words, any solution of class C 2 of the system J]} verifies the second order Euler- 
Lagrange equations (i.e., the jet geometric dynamics associated to the ODE 
system ([T])) produced by the jet least squares Lagrangian function ©: 



Remark 1 Conversely, the preceding statement is not true. In other words, 
there exist solutions of the Euler- Lagrange ODE system of second order ([3]), 
which are not global minimum points for the jet least squares energy action 
£ ODEs , that is which are not solutions for the initial jet dynamical system |T]). 

In such a context, we consider that we may regard the jet least squares La- 
grangian function JLS 00 ^ as a natural geometrical substitut on the 1-jet space 
J 1 (K, M. n ) for the dynamical system ((T|). But, an entire single-time Lagrange 
geometry on the 1-jet space J 1 (R, W l ) (in the sense of nonlinear connection, 
generalized Cartan linear connection, d-torsions, d-curvatures, jet "electromag- 
netic" d-field and jet "electromagnetic" Yang-Mills energy), geometry which is 
produced only by the jet least squares Lagrangian function JLS ODEs (via its 
Euler-Lagrange equations (J3J), is now completely done in the book [3J. For such 
a reason, we introduce the following concept: 

Definition 2 Any kind of geometrical object on J 1 (R, M. n ), which is produced by 
the jet least squares Lagrangian function JLS ODEs (via its second order Euler- 
Lagrange equations (J3J ) is called a geometrical object produced by the jet 
dynamical system ([!]). 

Let us consider the Jacobian matrix 



In such a context, we give the following geometrical result (which is proved in 
the book [3j): 

Theorem 3 (1) The canonical nonlinear connection on J 1 (M, M. n ) pro- 
duced by the jet dynamical system has the local components 




V i = 1, n. 



(3) 





where 




2^(2) 



T J(X)]. 
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(2) All adapted components of the canonical Cartan linear connection 
CpODEs p roc i uce( i fry th e jet dynamical system (JXJ are zero. 

(3) The torsion tensor T ODEs of the canonical Cartan linear connection 
cr ODEs produced by the jet dynamical system (JXJ) has the adapted 
components (k = l,n) 



(i)fc 



d a (1) 



f)2 jr(j) 
6 A (l) 

dx k dxi dx k dx % 



gN(i) 
dx k 



(4) All adapted components of the curvature tensor R ODEs f the canon- 
ical Cartan linear connection cr ODEs produced by the jet dynamical 
system ^ cancel. 

(5) The geometric "electromagnetic" distinguished 2-form produced 
by the jet dynamical system $Q has the expression 



where 

and we have 



5x\ = dx\ + N ( ( ^ )k dx k , V i = l,s 



dX 



(0 

(i) 



dX 



dxi 



dx l 



-N (1 



(!)• 



(6) The jet geometric "electromagnetic" Yang- Mills energy produced 
by the jet dynamical system ([1]) is given by the formula 



n— x n 

EYM ODEs (a;) = £ £ 



i=l j=i+l 



2 1 

= - • Trace 
2 



p(!) . t fW 



Remark 4 XTie adapted components °f ^e "electromagnetic" (d-) 2-form 

jrODEs p ro d uce d fry ffoe j e f dynamical system ([I]) verify the following geomet- 
rical Maxwell equations: 

X! F Siiifc =0 > 

where Y]u j k y represents a cyclic sum and we have 



dF 



(i) 



(i)j||fe g x k 

For more details in the jet "electromagnetic" topic, please consult the work [3J. 



4 



Remark 5 The jet geometric Yang-Mills energy EYM™ 5 coincides with the 
norm of the skew- symmetric "electromagnetic" matrix F^ 1 ' in the Lie algebra 

o(n) = L{0{n)) = { A e M n (K) | A + T A = 0} , 

where 0(n) = { v4 € M n (R) | ^4 • T A = /„} is ifte corresponding Lie group of 
the orthogonal matrices. 

3 Three-dimensional mathematical model used 
for studying the tibio-femoral dynamics 

The femur and tibia are modeled as two rigid bodies. Cartilage deformation is 
assumed relatively small compared to joint motions and not to affect relative 
motions and forces within the tibio-femoral joint. Furthermore, friction forces 
will be neglected because of the extremely low coefficients of friction of the 
articular surfaces. Hence, in this model, the resistance to motion is essentially 
due to the ligamentous structures and the contact forces. The menisci were not 
taken into consideration in the present model. For more biomechanical details, 
see Hefzy and Abdcl- Rahman [5]. 

The joint coordinate system, which was introduced by Grood and Suntay [5], 
is used to define the rotation and translation vectors that describe the three- 
dimensional patella-femoral and tibio-femoral motions. This joint coordinate 
system consists of an x-axis that is fixed on the femur (i is the unit vector 
directing the rc-axis), a z'-axis that is fixed on the tibia (k is the unit vector 
directing the z'-axis), and a floating axis perpendicular to these two fixed axes 
(i x k is the unit vector directing the floating axis). The Grood-Suntay's rotation 
vector includes three angular components: 

1. a is the flexion- extension angle that characterize the rotation occured 
around the femoral fixed axis Of x; 

2. /3 — 7r/2 is the adduction- abduction angle for the right knee, that charac- 
terizes the rotation arround the floating axis (this is the case studied by 
us). For the left knee, the adduction-abduction angle is given by it/2 — (3. 
In both cases the angle /3 is called the varus-valgus angle. 

3. 7 is the internal- external tibial angle that characterize the rotation which 
occurs about the tibial fixed axis Otz'. 

Using this joint coordinate system, the rotation vector (describing the ori- 
entation of the tibial coordinate system with respect to the femoral coordinate 
system) is defined as: 

8 — —ai — P \ i x k'\ — jk . 
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This rotation vector can be transformed to the femoral coordinate system, 
and then, it can be differentiated with respect to time to yield the angular 
velocity and angular acceleration vectors of the tibia with respect to the femur. 

In this analysis, it is assumed that the femur is fixed while the tibia is moving. 
The transformations of these two Grood-Suntay coordinate systems are given 
by (see 0,0) 

r = r OT + r' = t 0t + 9t-r t= o, 

where the vector fo T i s the position vector which locates the origin of the tibial 
coordinate system with respect to the femoral coordinate system (the tibial 
origin Ot is considered in the center of tibia) , ft=o describes the position vector 
of an arbitrary point P with respect to the tibial coordinate system (at the 
initial moment t — 0; note that the vector ft=o is constant in time), and r is 
the position vector (at an arbitrary moment t) of the same point P with respect 
to the tibial coordinate system. As usual in the mechanics of the rigid bodies, 
the position vector r' of the point P is expressed by the rotation formula 

r' = $H-r t=0 , 

where the matrix of rotation 

(i?n Ru R13 \ 
R21 R22 R23 (4) 
R31 R32 R33 j 

has the entries (see [5], [5]) 

i?n = sin /3 cos 7, R\2 = sin /3 sin 7, i?i3 = cos,9, 
R21 = — cos a sin 7— sin a cos f3 cos 7, 
R22 = cos a cos 7— sin a cos j3 sin 7, R23 = sin a sin /?, 
R31 = sin a sin 7— cos a cos /3 cos 7, 
R32 = — sin a cos 7— cos a cos j3 sin 7, R33 = cos a sin (3. 

The rotation matrix 9\ is generated by multiplying the three matrices generated 
by the rotations arround the axes, namely: 

1. a rotation of angle a arround the axis Of x (from the top of the vector i 
the rotation is clock- wise): 

/ 1 

Ko F i, a '■ = I cos a. sin a 
\ — sin a cos a 

7T 

2. a rotation of angle — — f3 arround the floating axis (from the top of the 
vector i x k the rotation is clock- wise): 

sin (3 cos /3 

n floating axis, /3 : = | 1 

- cos ft sin /3 
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3. a rotation of angle 7 arround the axis O^z' (from the top of the vector k 
the rotation is clock- wise): 

(cos 7 sin 7 
— sin 7 cos 7 
1 

In other words, we have the matrix equality 

— ^-Oi^x, a ' i** floating axis, /? ' I^Otz' . 7* 



The equations that govern the three-dimensional motion of the tibia with 
respect to the femur are the second order differential Newton's and Eulcr's 
equations of motion. 

The classical Newton's equations are written in scalar form, with respect to 
the femoral fixed system of axes, as: 

( F^ + G x =mi OF 
\ F^ + G y = myo F 
{ F^ + G z =mz OF , 

where 

• F ext — F° xt . i + Fy Xt ■ j + F° xt • k is the sum of all external forces applied 
to the tibia-femur (contact forces in knee, ligamentous forces in knee etc.); 

• G — G x ■ i + G y ■ j + G z ■ k is the weight of the leg of mass m; the mass 
of the leg was taken in this experiment as m — 4.0 kg (see Hefzy and 
Abdel- Rahman [B]). 

• o-Of = %o F ' * + Vo F ' j + zo F ' k is the acceleration of the center of mass of 
the leg (which is fixed in the origin of the femoral coordinate system Of)- 

From the perspective of the rigid body mechanics, the Euler's equations of 
motion are written in the scalar form as: 

I x > x >(b x > + (I z > zl - I v /y>)u v 'UJ z > = M^f* 

< Iy'y'LUy, + {I x , x . - I Z > Z >)U Z >U X > = (5) 
^ I Z > Z 'U Z > + {Iy>y> ~ I X ' X >)U X 'U)y> = , 

where 

• 57 = uj x i ■ i + ujyi ■ j' + uj z i ■ k is the angular velocity vector of the tibia 
with respect to the femur; 

*_ ._/._/.—/ 

• ui = ui x > ■ i + LUyi ■ j + uj z i ■ k is the angular acceleration vector of the tibia 

with respect to the femur; 
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• M cxt = MJ?* • t + My?* 1 ■ j' + M°, xt • k' is the sum of the moments (the 
torques) of all external forces acting on the tibia around the x'-, y'-, and 
z'-axes; 

• Ix'x'-, ly'y' and I z i z i are the principal moments of inertia of the leg about 
its centroidal principal system of axes. The inertial parameters were esti- 
mated using anthropometric data as 

I x , x , = 0.0672kg- m 2 , = 0.0672kg- m 2 , I z , z , = 0.005334kg- m 2 . 

Note that, in this analysis, the leg was assumed to be a right cylinder. 
For more details, please see Hefzy and Abdel- Rahman J5], and references 
therein. 

Remark 6 If we use the matrix equality 

( t, k') = (l 3, *)•*, 

then we deduce that the rotation vector 9 takes in the femoral system of axes the 
following form: 

6 = (— a — 7cos/3) • i + (— f3 cos a — 7 sin a sin /3) • j + (6) 
+ (J3 sin a — 7 cos a sin /3) ■ k. 

Therefore, by differentiating the rotation vector 6, we find the angular velocity 
vector _ 

d0 — — — 

U> = — = U) x ■ 1 + LJ y ■ J + U) z ■ K. 

Now, by rotating the angular velocity vector ZJ, via the formula 

( u x /, u v >, oj z > ) = ( w x , u y , u> z ) ■ £H, 

we obtain the equalities (see also [6]^: 

w x i — — d sin /3 cos 7 — d/3cos/3cos7 + 07 sin f3 sin 7 + fj sin 7 + /37COS7, 
LJ y / = — a sin fj sin 7 — a/3 cos (3 sin 7 — cey sin f3 cos 7 — f3 cos 7 + /?7 sin 7, 
w z ' = —a cos /3 + a/3 sin /3 — 7. 

4 Approximate values of the components of the 
total external torque and angular velocity 

Using the three-dimensional mathematical model for the tibio-femoral motion 
due to Luh et al. (see [S]), Apkarian et al. (see pQ) practically measured in a lab 
frame (using some specialized devices) the total external torques acting on the 
tibia of a subject (A) having the stride length 1.41 m and the speed of walking 
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1.21 m/s. Obviously, the time of the gait cycle for the subject (A), studied in 
the paper pQ, is T = 1.1652 s. 

Let us consider for the subject (A) the following gait cycle intermediate 
moments: 

t = 0(0%), h = T/4 = 0.2913 (25%), t 2 = T/2 = 0.5826 (50%), 

t a = 3T/4 = 0.8739 (75%), t 4 = T= 1.1652 (100%). 

Looking now at the three graphs of the knee torque components (computed 
in N • m in the paper pQ, pp. 153), we see that the total external torque vector 
of the knee 

M cxt = M£ xt • i + M* xt • J + M| xt ■ k := 

-^valgus ' i -^extension ' j H~ -^internal ' K 

has the following approximate interpolation components: 



M^(t ) w 7.5, 


A/, cxt (t ) « 7.5, 


M™(t )=Q, 


M^(h) ps -40, 


M y ext (*i)=0, 


M 2 oxt (*i) « 5, 


M x ext (* 2 ) « -15, 


M* xt (t 2 ) = 0, 


M z oxt (t 2 ) =0, 


M^(t 3 ) = 0, 


M£ xt (t 3 )=0, 


Mr*(t 3 ) « -5, 


M£ xt (t 4 ) « 7.5, 


M y ext (t 4 ) « 15, 


M 2 oxt (t 4 ) = 0. 



At the same time, looking at the approximate values of the three knee angles 
appearing in the graphs from [1 (pp. 150), we observe that the rotation vector 

= 6 x ■ l -\- 6 y ■ j -\- 9 z ■ k '.= $ varus ' i $ flexion ' j ^external ' k 

has the interpolation components (the angles are presented here in radians; from 
the top of the vectors i, j, fe, these angles are measured counter-clock- wise): 



X (U>) = (= 0°), 


9 v (t ) ps -0.0872 (= -5°), 


6 x (ti) ps -0.0872 (= -5°), 


9 y (ti) ps -0.3490 (= -20°), 


9 x {t 2 ) ps -0.0872 (= -5°), 


y (t 2 ) ~ -0.3490 (= -20°), 


x {h) ps -0.1745 (= -10°), 


9 y (t 3 ) ps -1.1344 (= -65°), 


e x {u) = o (= o°), 


9 y (U) ps -0.0872 (= -5°), 



9 z (h) R 


i -0.1745 (= 


-10°), 




i -0.0872 (= 


-5°), 




i -0.1745 (= 


-10°), 


0»(* 3 )R 


i -0.0872 (= 


-5°), 




i -0.1745 (= 


-10°). 
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Remark 7 Solving numerically the angle system ( see the rotation vector ^ ) 

!-a - 7 cos/? = 6 X 
—j3 cos a — 7 sin a sin (3 = 9 y 
(3 sin a — 7 cos a sin (3 — 9 Zl 
we find the following intermediate Grood-Suntay's angles for knee: 



a(t ) « -0.5786, 


/3(io) « 0.1684, 


7(t ) w 0.5869, 


a(ti) « -0.0760, 


/3(h) w 0.3546, 


7(ii) « 0.1740, 


a(t 2 ) ~ -0.1851, P{h) w 0.3751, 


7(ia) « 0.2927, 


a(t 3 ) ~ 0.0859, /3(t 3 ) « 1.1227, 


7 (f 3 ) « 0.2044, 


a(t 4 ) « -0.5786, /3(i 4 ) « 0.1684, 


7(i 4 ) « 0.5869. 



Now, using the Lagrange polynomial of interpolation generally described by 



yo • 


x — 


Xi 


X — 


X2 


X — 


X3 


x — 


£4 


Xq - 


Xl 


Xq - 


x 2 


Xq - 


X3 


Xq ~ 


— + 

X4 


+Vi 


X 


- x 


X 


- x 2 


X 


- ^3 


X 


— X4 


X\ 


- Xq 


Xl 


- x 2 


Xl 


- x 3 


Xl 


— x 4 


+ V2 


X 


- Xq 


X 


- Xl 


X 


- X3 


X 


— X4 


X2 


- Xq 


X2 


- Xl 


X2 


- x 3 


X2 


— x 4 


+ V3 


X 


- Xq 


X 


- Xl 


X 


- x 2 


X 


— X4 


X3 


- Xq 


X3 


- Xl 


X3 


- X2 


X3 


— x 4 


+Vi 


X 
X4 


- Xq 

- Xq 


X 

X4 


- Xl 

- Xl 


X 
X4 


- x 2 

- x 2 


X 
Xi 


- X3 

- x 3 



we can construct the approximate Apkarian's angle functions x (t), 9 v (t) and 
9 z (t). Afterward, by differentiating the angle functions 9 x (t), 9 y (t) and 9 z (t), we 
get the components of the angular velocity in our gait cycle moments (rad/ s): 



uj x (t ) w -1.0981, 


u y (t ) « -5.6920, 


w z (i ) « 1.5984, 


u> x (h) w 0.0999, 


u>y(h) « 1.1983, 


« -0.3997, 


u x (te) ~ -0.1998, 


w„(ta) « -1.7974, 


w,(ia) « -0.0002, 


wx(i 3 ) « -0.1997, 


w„(i 3 ) w -2.0970, 


w,(i 3 ) « 0.3993, 


u x {U) rj 1.8975, 


Wj,(i 4 ) « 12.8820, 


w*(t4) « -1.5983. 



Consequently, via the rotation formulas 

X = i?nX + i?2lY + i?3iZ, 

Y = i? 12 X + i? 22 Y + R 32 Z, 

Z = i2i 3 X + i?23Y + -R33Z, 

where R\j are given by the expressions from l|4"p. we find the following interme- 
diate vector components (with respect to tibia coordinate system): 
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■ — components of the total external torque: 



Mg*(to) w 0.9361, 


M y ext (to) ps 8.1637, 


M°r(t ) ps 6.7065, 


M^*(ti) ps -18.3490, 


M^*(ti) ps -2.8400, 


M z c , xt (t!) ps -35.7800, 


M| xt (t 2 ) w -5.2618, 


M* xt (t 2 ) « -1.5857, 


M^(t 2 ) ps -13.9570, 


M* xt (t 3 ) » 2.0263, 


M^(t 3 ) ps 0.8581, 


M z c , xt (t 3 ) ps -4.4898, 


M| xt (t 4 ) ps 0.8255, 


M£ xt (i 4 ) « 15.6310, 


M£ xt (t 4 ) ps 6.0191. 


■ — components of the angular velocity: 




uvOo) « -1-6519, 


uv(to) ps -5.7721, 


w*'(*o) « -0.3366, 




c«v(ti) w 0.2847, 


w„'(ti) ~ 1.2324, 


ps -0.0762, 




uv(i 2 ) w 0.1451, 


Lo y (h) ps -1.8010, 


w«'(*2) ~ -0-0647, 




w-e/fts) « 0.1623, 


aty(t 3 ) ps -2.1350, 


w^(t 3 ) « 0.1098, 






uv(i 4 ) « 1.6574, 


wj,/(t 4 ) ps 13.0050, 


uv(t 4 ) « 0.4656. 





Now, using the preceding Tables and the well-known method of multiple re- 
gression from Statistics, by numerical computations, we get the following linear 
approximations for the torque components M™', M™' and M| xt , with respect 
to the angular velocity components u) x i, uj y > and u z > (via the classical method 
of least squares) : 



16.9430 • u x , - 0.5003 • uy + 80.9290 • u x > - 3.0709, 
15.1720 • LJ x f + 1.5740 • uy + 34.8270 • u z , + 3.7511, (7) 
39.7610 • u x > + 0.9071 • uy + 140.8000 • uy - 7.1266. 

5 From Euler's equations of motion of knee to 
jet "electromagnetic" Yang-Mills energy 

Note that, in our jet geometrical approach, we can regard the Euler's equa- 
tions of motion ((SJ as a dynamical system on the 1-jet space J 1 (R, M 3 ). The 
coordinates on the 1-jet space J 1 (R, M 3 ) are considered as being 

In such a context, taking the particular case from the Hefzy - Abdel-Rahman's 
paper [6], when we have 

W = Iyy = 0.0672, I z , z , = 0.0053, 

and setting the components of the total external torque by the expressions (|T|), 
then the Euler's equations of motion (|SJ) for the subject (A) can be rewritten in 



M ext ps - 
M cxt ps - 
M c , xt ps - 
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the following particular form: 
civ = 0.9211 • LU qi 



- 252.1279 • uv - 7.4449 • ay + 



(8) 



+1204.3005- uv - 45.6979 

u) y , = -0.9211 • uj x , -uj z , - 225.7738 • u x > + 23.4226 -uy + 

+518.2589 -uj z , +55.8199 

Cj z , = -7502.0754- u x , + 171.1509 • ay+ 

+26566.0377 • u z , - 1344.6415. 

Consequently, via the Theorem [3] applied to the first order ODE system ©, 
we assert that the Lagrangian geometrical behavior on the 1-jet space J X (R, R 3 ) 
of the above knee dynamical system can be described by the following result: 



Corollary 8 (1) The canonical nonlinear connection on J 1 (M,R 3 ) pro- 
duced by the jet dynamical system of knee ([SJ has the local compo- 
nents 

where the matrix of the spatial nonlinear connection 



( N (1) /V (1) 

' iV (l)l JV (1)2 



N(i) = 



/v (2) 

^(1)2 
r( 3 ) 



(3) 



^(1)2 JV(1)3 / 



/ 

* (1)2 
\ *(1)3 



-*(1)2 



(2)3 



"*(1)3 ^ 
* (2)3 

o J 



has the entries 

0.9211 • u z , + 109.1644 



* (1)2 



F (! 



(1) 
)3 



= 0.4605 -a; 



4353.1879, 



(2)3 



0.4605 -av + 173.5540. 



(2) All adapted components of the canonical Cartan linear connection 
CI 1 produced by the jet dynamical system of knee (JSj) are zero. 

(3) The torsion tensor T o/ </ie canonical Cartan linear connection CT pro- 
duced by the jet dynamical system of knee has as adapted com- 
ponents the entries of the following torsion matrices: 



dN (1 



(i)i 



(i) 



T(i 



(1)2 



gN(i) 

duly' 



( 



V o 
/ 




0.4605 
-0.4605 



-0.4605 \ 

\ 0.4605 j 
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dN (1 



(1)3 



(1) 



3uj z i 



-0.9211 \ 
0.9211 
j 



(4) All adapted components of the curvature tensor R of the canonical Car- 
tan linear connection CT produced by the jet dynamical system of 
knee ([8]) cancel. 

(5) The geometric "electromagnetic" adapted components produced 
by the jet dynamical system of knee © are the entries of the matrix 



Hi) 



-N 



(i) 



V 



o 

pw 

* (1)3 



«,(!) 
(1)2 



F (1) 

* (2)3 



*(1)3 ^ 
* (2)3 



(6) The jet geometric "electromagnetic" Yang-Mills energy produced 
by the jet dynamical system of knee (jHJ) is given by the formula 



EYM kncc K,,av, Wz ,) 



(1)2 



+ 



(1) 



(2)3 



Proof. The Euler's equations of motion ([5]) represent a particular case of the 
jet hrst order ODE system ([T]) for n = 3 and X = (xl\(uj X ',ujy>,u) z >)) , 

V / i=l,3 



where 



X^iujx^ujy^ujz,) = 0.9211 • u>y, -uj z , - 252.1279- u> x > - 7.4449 • u y > + 
+1204.3005 • u~> - 45.6979, 



-0.9211 ■ uj x , ■ uj z , - 225.7738 • u x , + 23.4226 • w„ 



-518.2589 -w 2 



55.8199, 



,(3) 



InK,^,^') = -7502.0754. 



171.1509 -a; 



+26566.0377 ■ ui z , - 1344.6415. 



Consequently, using the formulas 



b (1)2 



duj x i 



dojy' 



p(V 

' '(1)3 



dX {1) dX {1) 



Ouj x > 



(2)3 



9w, 



9w u 



F« = -N (1 



(i)' 



some direct computations lead us to what we were looking for. 
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6 Yang-Mills energetic surfaces of constant level 
produced by the jet dynamical system of knee 

Using the notations ui x i := X, u> y > :— Y, u z > := Z, then the jet geometric 
" electromagnetic" Yang-Mills energy produced by the dynamical system of knee 
© (associated to the subject (A)) takes the form 

EYM kncc (I, Y, Z) = [0.92112 + 109.1644] 2 + [0.4605Y + 4353. 1879] 2 + 

+ [-0.4605X + 173.5540] 2 . 

Now, let us consider the jet "electromagnetic" Yang-Mills energetic surfaces 
of constant level, which are produced by the dynamical system of the knee of 
the subject (A). These energetic surfaces are defined by the implicit equations 

S kncc . E YM knoc (A, Y, Z) = k, 

where k > is a given constant real number. By direct numerical computations, 
we obtain the following approximate implicit equations: 

S^ nco : 0.2120 {X'f + 0.2120 (Y'f + 0.8484 (Z'f =k^ 

E knee . ( X 'f , ( Y 'f , ( Z 'f _ 1 = Q t. U Q 

k ' 4.7169/c 4.7169fc 1.1786fc ' r ' 

where we have done the spatial translation 

X' = X- 376.8816, Y' = Y + 9453.1767, Z' = Z + 118.5152. 

In conclusion, it follows that the jet "electromagnetic" Yang-Mills energetic 
surface of constant level k > is: 

1. for k = 0, it is the point C(376.8816, -9453.1767, -118.5152). 

2. for k > 0, it is an ellipsoid centered in the point C, whose semi-axes are 
parallel with the axes OX, OY and OZ, and they have the lengths 

a = b = \/4.7169fc, c = Vl.l786fc. 

Open problem. Find the biomechanical interpretations for the ellipsoid 
shapes (oblate spheroids) of the jet Yang-Mills energetic surfaces of constant 
level produced by the dynamical system of knee ([5]) . 
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